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Abstract 

Push-pull networks are ubiquitous in signal transduction pathways in both prokaryotic and 
eukaryotic cells. They allow cells to strongly amplify signals via the mechanism of zero-order 
ultrasensitivity. In a push-pull network, two antagonistic enzymes control the activity of a protein 
by covalent modification. These enzymes are often uniformly distributed in the cytoplasm. They 
can, however, also be colocalized in space, for instance, near the pole of the cell. Moreover, it is 
increasingly recognized that these enzymes can also be spatially separated, leading to gradients 
of the active form of the messenger protein. Here, we investigate the consequences of the spatial 
distributions of the enzymes for the amplification properties of push-pull networks. Our calculations 
reveal that enzyme localization by itself can have a dramatic effect on the gain. The gain is 
maximized when the two enzymes are either uniformly distributed or colocalized in one region 
in the cell. Depending on the diffusion constants, however, the sharpness of the response can be 
strongly reduced when the enzymes are spatially separated. We discuss how our predictions could 
be tested experimentally. 
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Synopsis 

Living cells continually have to respond to a changing environment. To this end, they 
do not only have to detect environmental signals, but also amplify them. In living cells, 
signals are often amplified in so-called push-pull networks. In a push-pull network, two 
enzymes control the activity of a protein in an antagonistic manner. A well-known example 
is a network in which a kinase phosphorylates a messenger protein, while a phosphatase 
dephosphorylates the same protein. While it has long been assumed that the enzymes 
are uniformly distributed in the cytoplasm, it is increasingly becoming clear that in many 
systems one or both of the enzymes are localized in space, for instance near the cell pole. 
If the enzymes are spatially separated, then spatial gradients of the messenger protein can 
form, and recently a number of these protein gradients have been observed experimentally. 
We study by numerical calculations how the amplification properties of push-pull networks 
depend upon the spatial distribution of the enzymes. We find that the gain is maximized 
when the enzymes are either uniformly distributed or colocalized in space. Depending upon 
the diffusion constants, however, the sharpness of the response can be strongly reduced when 
the enzymes are spatially separated. 
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Introduction 



Living cells are information processing machines. In order to process information reliably, 
signals often need to be amplified. To this end, cells can employ a variety of amplification 
mechanisms. Signals can be amplified via positive feedback, cooperative binding of signaling 
molecules to receptors, or interactions between receptor molecules jlL Another principal 
mechanism for signal amplification is zero-order ultrasensitivity j^, 3|. This mechanism 
operates in so-called push-pull networks, which are omnipresent in both prokaryotes and 
eukaryotes. In a push-pull network, two enzymes covalently modify a component in an 
antagonistic manner (see Fig. [1]). One well-known example is a network in which a kinase 
phosphorylates a component, and a phosphatase dephosphorylates the same component. If 
both enzymes operate near saturation, then the modification reactions become zero order, 
which means that the reaction rates become insensitive to the substrate concentrations. 
Under these conditions, a small change in the concentration of one of the two enzymes (the 
input signal), will lead to a large change in the concentration of the modified protein (the 



output signal 



in detail Q, li 



2,3]. The amplification properties of push-pull networks have been analyzed 



In these studies, however, it is assumed that the antagonistic 



enzymes are uniformly distributed in space. Yet, it is increasingly recognized that in many 
systems one or both of the two antagonistic enzymes are localized in space, for instance at 
the cell pole. Here, we address the question how the spatial distribution of the antagonistic 
enzymes affects the amplification properties of push-pull networks. 

If the two antagonistic enzymes are separated in space, then gradients of the messenger 

Il2 . [13I ] . Recently, a number of protein gradients have been ob- 



protein can form 
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served experimentally in both prokaryotic and eukaryotic cells. For example, in Escherichia 
coli cells, the kinase CheA and the phosphatase CheZ control the phosphorylation level of 
the messenger CheY, which transmits the chemotactic signal from the receptor cluster to 
the flagellar motors. In wild-type cells, the kinase and the phosphatase are both localized 
at the receptor cluster jl^j ]. and, as a result, the steady-state concentration profile of CheY 
is uniform [h|. However, in E. coli mutants, where the phosphatase is distributed in the 



cytoplasm, gradients of CheY have recently been observed 



Other examples of protein 



gradients include Caulobacter, in which MipZ gradients guide chromosome segregation and 



division site selection 



151 ] . In eukaryotic cells, gradients of Ran, Stathmin, and HURP pro- 
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teins aid in the formation of the mitotic spindle by providing directional cues for microtubule 
growth If], 17, IsiQl- Moreover, in eukaryotic cells, the kinases in the mitogen-activated 
protein kinase (MAPK) cascade often bind to scaffold proteins, while the phosphatases are 



distributed in the cytoplasm 



20] . This will lead to concentration gradients of the acti- 



vated kinases, which can become particularly important if the scaffolds are located near the 
membrane. 

In this study, we compare the amplification properties of a canonical push-pull network, 
where all components are uniformly distributed in space, with those of a network in which 
the enzyme that provides the input signal is localized at one end of the cell, while all the 
other components can freely diffuse through the cell. In the latter case, the concentration 
profile of the messenger — the output signal — is non-uniform. Previous studies have focused 
on the time-dependent concentration profiles of the messenger , ll3| and on the 'control' 



of diffusion over protein fluxes 21| in similar systems. Here, we examine the effect of the 
spatial distribution of the enzymes on the amplification properties of push-pull networks. 
To this end, we compute for both systems the steady-state input-output relations. Our 
analysis reveals that the spatial distribution of the enzymes can have a dramatic effect on the 
capacity of push-pull networks to amplify input signals: the maximum gain of the network 
in which one enzyme is localized at one end of the cell, while the other is not, can be much 
lower than that of the network in which the components are uniformly distributed in space. 
Importantly, this effect occurs over a range of diffusion constants, protein concentrations, 
and enzymatic activities that is typical for living cells. 

In the next section, we introduce the push-pull network. In the Results section, we first 
present the input-output relations for both networks. We show that the gain can be much 
reduced when the enzymes are spatially separated, and demonstrate that the magnitude of 
this effect depends upon the diffusion constants of the diffusing components. To elucidate 
the dose-response curves, we discuss in the subsequent sections the spatial concentration 
profiles in both the low and high activation limits. This analysis reveals that the maximum 
gain in the non-uniform system is reduced, because the response of the network depends 
on the position in the cell. Interestingly, the calculations also show that separating the 
enzymes in space does not only attenuate strong signals by limiting the maximum response, 
but can also enhance the propagation of weak signals. 



4 



Methods 

The push-pull network 

A push-pull network consists of two Michaelis-Menten reactions (see also Fig. [1]): 

E a + X^E a X^E a + X* (1) 

E d + X* % E d X* ^ E d + X (2) 

fcs 

Here, E a is the activating enzyme that provides the input signal, and Ed is the deactivating 
enzyme. The substrate X is the unmodified messenger that serves as the detection compo- 
nent and X* is the modified messenger that provides the output signal; E a X denotes the 
activating enzyme bound to its substrate X and EdX* is the deactivating enzyme bound to 
its substrate X*. 

If all the components are uniformly distributed in space, then the chemical rate equations 
that correspond to this network are: 

d[X*\ 



dt 

at 

dt 
<9[E a X] 

dt 

d[E d ] 



dt 

<9[E d X* 



fc 3 [E a X] - k 4 [E d }[X*] + k 5 [EdX*} (3) 

^[EdX^-A^EjKR^EaX] (4) 

(fc 2 + fc 3 )[E a X]-A; 1 [E a ][X] (5) 

k 1 [E a ][X]-(k 2 + h)[E a X] (6) 

(fc 5 + A; 6 )[EdX*]-A; 4 [Ed][X*] (7) 

h[E d ][X*}-(k 5 + k 6 )[E d X*} (8) 



dt 

Here, [...] denotes the concentrations of the species. The steady-state input-output curve of 
this network can be obtained analytically 

We will compare the behavior of this network to that of a network in which the activating 
enzyme E a is located at one pole of the cell, while the other components can freely diffuse 
in the cytoplasm. The cell is assumed to be cylindrically symmetric. Since we are interested 
in the mean concentration profiles, it is meaningful to integrate out the lateral dimensions y 
and z. We thus consider a simplified ID model, with concentrations as a function of x only. 
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This leads to the following reaction-diffusion equations: 

3£! = + k 3 [E a X]5(x) - k 4 [E d ][X*] + 

k 5 [EdX*] (9) 
ffl = D^H + A; 6 [E d X*]-A; 1 [E a ][X]5(x) + 

k 2 [EaX]5(x) (10) 



g[gj 

at 

<9[E a X] 
dt 

0[Ed] n <9 2 [E 



(fe + ^CEaXJ-fcitEJpCKO) (11) 
fc![E a ][X](0)-(A; 2 + A;3)[E a X] (12) 



jD ^fT ii + ( fc 5 + ^)[E d X*]- 



<9t <9x 2 

fc 4 [E d ][X*] (13) 
<9[E d X*] _ ^^[EdX* 
dt dx 2 

(£; 5 + fc 6 )[E d X*] (14) 



The components E a and E a X are localized in the membrane at one end of the cell; the unit 
of their concentrations is number of molecules per area. The other components diffuse in the 
cell. Their concentrations, which are in units of number of molecules per volume, depend on 
the position x in the cell, where x measures the distance from the pole at which E a and E a X 
are localized; only in Eqs. [TT] and [12] is the x dependence explicitly indicated to emphasize 
that the E a -X association rate depends on the concentration of X at contact. Zero- flux 
boundary conditions are imposed at both cell ends. The steady-state input-output relations 
of the network described by Eqs. IMT41 were obtained numerically by discretizing the system 
on a (ID) grid and propagating Eqs. [9lfT4lin space and time until steady-state was reached. 
We consider a cell with the typical dimensions of an E. coli cell: the length of the cell, 



L, is thus on the order of 3/im 10] . We assume the same diffusion constants for all the 
components that can diffuse in the cytoplasm. This is for reasons of simplicity; it is not 
essential for the main conclusions of our work. To focus on the effect of enzyme localization 
on the input-output relation, we assume for both networks that k\ = k±, k-i = k$, ks = kg] 
the Michaelis-Menten constants for the modification and demodification reactions are thus 
the same: Km,s, = (&2 + &3V&1 = ^M,d = (k 5 + k 6 )/k4. To compare the two networks on equal 
footing, the total concentration of activating enzyme, [E a ]x = [E a ] + [E a X], was chosen such 
that [E a ]^ u = L[E a ]^, where [E a ]^ u is the concentration (per unit area) in the non- uniform 
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system and [E a ]^ is the concentration (per unit volume) in the spatially uniform network. 
This choice ensures that the total number of activating enzyme molecules in the whole cell 
is the same for both systems. In what follows, we will report [EJt = [E a ]x- 

In the calculations, we vary the concentration of the activating enzyme, E a , which 
is the input signal. The total concentration of the deactivating enzyme, Ed, is kept 
constant at [E^t = 0.5/iM; the rate constants are fixed at k\ = k<± = 10 8 M~ 1 s _1 , 
^2 = ^5 = 25 s -1 , ks = kg = 25 s -1 , corresponding to Michaelis-Menten constants of 
K M = Kyi A = i^M,d = 0.5/iM. We will study extensively the effect of changing the 
diffusion constant D and the total substrate concentration [S]t = [X]t + [X*]t, where 
[X] T = [X] + [E a X]/L is the total concentration of X and [X*] T = [X*] + [E d X*] is the 
total concentration of X*. Their base-line parameters, however, are: D = 10/im 2 s -1 and 



[S]t = 20//M. The magnitude of the diffusion constant 



as well as the values of the 



Michaelis-Menten constants, enzyme concentrations, and substrate concentrations, are 



typical for prokaryotic |23j] and eukaryotic cells [4|. 



Results 

The input-output relation 

Goldbeter and Koshland showed that if the antagonistic enzymes in a push-pull net- 
work operate near saturation (see Fig. [1]), a small change in the concentration of the 
activating enzyme E a can lead to a large change in the output, the modified messenger 
X* Q]. The enzymes become more saturated with substrate when either the Michaelis- 
Menten constants Ku,a, and Ku,d decrease, or the total substrate concentration [S]t = 
[X] + [E a X]/L+ [EjX*] + [X*] increases. Fig. [2^ shows the steady-state input-output relation 
for a push-pull network in which all the components are uniformly distributed in space, for 
different substrate concentrations. It is seen that as the substrate concentration is increased, 
the sharpness of the response is drastically enhanced. This is the hallmark of the mechanism 
of zero-order ultrasensitivity. 



In many systems, such as the bacterial chemotaxis network of E. coli 10j, the two antag- 
onistic enzymes are colocalized at the same pole, while the detection component X and the 
messenger X* can diffuse through the cytoplasm. While the time-dependent response curves 
of such a network will differ from those of the two networks considered here, the steady-state 
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dose-response curves will be identical to those of a network in which all the components are 
homogeneously distributed in the cytoplasm. The response curves shown in Fig. [2^, thus 
also pertain to push-pull networks in which the two enzymes are colocalized at one end of 
the cell, while their substrates freely diffuse in the cytoplasm. Indeed, also in these networks 
the mechanism of zero-order ultrasensitivity can strongly amplify input signals. 

Spatially separating the enzymes reduces the gain. Fig [2b shows the dose-response curves 
for a push-pull network in which the activating enzyme E a is localized at one pole of the 
cell, while the other components diffuse in the cytoplasm. Three points are worthy of note. 
The first is that the maximum output signal, the concentration of the messenger X*, is 
much lower than that of the corresponding network in which all components are uniformly 
distributed in space (see Fig. [2]). In fact, while in the spatially uniform network, the fraction 
of modified substrate, [X*]t/[S]t, always approaches unity if [E a ]x/[Ed]x becomes large, in 
the non-uniform network the fraction of modified substrate saturates to a lower level: even 
when the concentration of activating enzyme is much higher than that of the deactivating 
enzyme, not all substrate X is converted into X*. The second point to note is that as the 
total substrate concentration decreases, the inflection point of the dose-response curve shifts 
to lower values of [E a ]x/[Ed]x- The last, and perhaps most important, point to note is that 
the sharpness of the response of the network is much weaker than that of the network in 
which the enzymes are either colocalized or uniformly distributed in space. The insets of 
Fig. [2] shows the logarithmic gain, g = din [X*]x/<91n [E function of [E a ]x/[Ed]x for 

both networks ([E<j]t is kept constant). It is seen that for both low [E a ]x/[Ed]x an d high 
[E a ]x/[Ed]x the gain is small and fairly similar for both networks, while for the symmetric 
networks considered here, at [E a ]x ~ [Ed]x the gain is maximal, but smaller for the network in 
which the enzymes are spatially separated. Hence, spatially separating the two antagonistic 
enzymes reduces the maximum gain of a push-pull network. 

The dose- response curves strongly depend on the diffusion constants. The extent to which 
the spatial separation of the opposing enzymes can change the response of the network 
depends on the magnitude of the diffusion constant of the components. This is illustrated 
in Fig. [31 This figure shows the input-output relation of a push-pull network where the 
activating enzyme is located at one end of the cell, while the other components diffuse freely 
in the cytoplasm, for different values of the diffusion constant. This network is in the zero- 
order regime: the total substrate concentration is large compared to the concentrations of the 



8 



enzymes and the Michaelis-Menten constants. Yet, for low values of the diffusion constants, 
the response is rather weak. As the diffusion constant increases, however, the sharpness of 
the response markedly increases. For D — > oo, the input-output relation approaches that of 
a push-pull network in which all components are either uniformly distributed in space, or 
colocalized in one region of the cell. 

Spatially separating the enzymes attenuates the propagation of strong signals, but can 
enhance the transmission of weak signals. Fig. [3] shows that in a zero-order network in 
which only the activating enzyme is localized at one pole of the cell, the concentration of 
X* decreases with decreasing diffusion constant when [EJt > [EJt, but increases with 
decreasing diffusion constant when [EJt < [EJt- This means that when the input signal is 
strong (high kinase activity), the response of a network in which the enzymes are spatially 
separated is weaker than that of a network in which the enzymes are either uniformly 
distributed or colocalized in space; conversely, when the input signal is weak (low kinase 
activity), the spatially non- uniform network can respond more strongly than a uniform 
network. Spatially separating the antagonistic enzymes will thus attenuate strong input 
signals, but can also amplify weak input signals. 

Mechanism: concentration gradients 

To explain the effect of enzyme localization on the amplification properties of push-pull 
networks, it is instructive to consider the effect of diffusion on the input-output relation: 
in the limit that D — > oo, the response of the network in which the activating enzyme is 
located at the pole, while the other is distributed in the cytoplasm, approaches that of a 
network in which the enzymes are either uniformly distributed in space or colocalized at the 
pole. The effect of diffusion on the response curves can be understood by considering the 
effects of diffusion on the spatial concentration profiles. 

In a push-pull network where the antagonistic enzymes are either uniformly distributed 
or colocalized in space, the steady-state spatial concentration profiles of the freely diffusing 
components are uniform across the cell. In a push-pull network where the two antagonistic 
enzymes are spatially separated, concentration gradients of the freely diffusing components 
can form. Fig. H] shows for a zero-order network in which the activating enzyme is located 
at one pole of the cell, while the other is not, the concentration profiles of X* and EdX*, for 
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three different (total) concentrations of the activating enzyme E a , [EJt- The concentrations 
of X* and E d X* are highest near the pole where X is activated, and decay in the cytoplasm 
where X* is deactivated. Moreover, the concentration profiles increase as [E a ]x increases. 
These gradients impose fundamental limits on the maximum gain of the system. 

To clarify the effect of diffusion on the concentration profiles and the input-output re- 
lations, it is instructive to recall that, in general, the spatio-temporal evolution of [X*] is 
given by the interplay of activation, deactivation and diffusion of X*: 

Here, J denotes the influx of X* into the system, while 7 denotes the deactivation rate of 
X* at position x. If the formation of the enzyme-substrate complexes is fast (see Text SI), 
then J and 7 are given by 

**> = ^ jiiw <i7) 

Here, [E a ] T = [E a ] + [E a X] and [E d ] T = [E d ] + [E d X*] are the total concentrations of E a 
and E d , respectively. Combining Eq. [13] with Eq. [14] reveals that the total concentration 
profile of Ed, [E^tO^)) is constant in space if, as assumed here, the diffusion constants of 
the enzyme Ed, and that of the enzyme bound to its substrate, EdX*, are the same. The 
synthesis rate of X* depends upon the concentration of X at contact and hence upon the 
concentration of X* at contact; similarly, the deactivation rate of X* at position x depends 
upon the concentration of X* at x. This is important to note, because, as we discuss 
below, the dose-response curves are determined by the sensitivities of the influx J and the 
deactivation rate 7 to changes in the substrate concentration. We will now first discuss 
the input-output relations of zero-order push-pull networks, and then briefly the response 
curves of push-pull networks that are in the linear regime. 

Push-pull networks in the zero-order regime 

Figs. HUH] show the concentration profiles and dose-response curves of push-pull networks 
that are in the zero-order regime. We now discuss the limits of weak and strong activation 
separately. 
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Weak activation 

We first consider the regime in which the concentration of the activating enzyme is lower than 
that of the deactivating enzyme, corresponding to Figs.[5|i-c. In the limit that [E a ]x *C [Ed]T, 
[X] will be large and [X*] will be small. As a consequence, E a is saturated with its substrate 
X, while Ed is not saturated with its substrate X*. Because E a is saturated, the influx of 
X* into the system is constant (i.e. independent of [X] and [X*]) and given by J = /^[EjxX 
(see Eq. [TBI) . Because E^ is unsaturated, the deactivation rate 7 is proportional to [X*]: 
7(x) = /i[X*](x), with \i = /c 6 /i~YM,d[Ed]T ( see Eq. ITT1) . This means that in this regime the 
deactivation rate per particle is constant. 

With the influx J being constant and the deactivation rate 7 being proportional to [X*] , 
Eq. [15] can be solved analytically (see Text SI) . Defining the characteristic decay length of 
X* to be A = y/D/ fi, then, if L ^> A, the solution is 

\K*](x) = ^exp(-x/\). (18) 

Eq. [26] reveals that when D increases, the profile decays more slowly, and the concen- 
tration of X* at contact decreases. When D increases, the X* molecules diffuse away from 
the pole more rapidly. Because in the regime considered here, namely [E a ] T <C [Ed]T, the 
influx of X* is constant and independent of D, the concentration of X* close to the pole will 
decrease when the molecules diffuse away more rapidly, while the concentration further away 
will increase. In fact, in this limit the total concentration of X*, [X*] cc11 , is independent of 
the diffusion constant; this can be verified by integrating Eq. [261 over the whole cell, which 
yields [X*] cc11 = J//i. The fact that the total concentration of [X*] is independent of the 
diffusion constant, means that the response of the network does not depend upon the spatial 
distribution of the enzymes. 

When [EJt increases, [X*] increases and [X] decreases. As a result, E a becomes less 
saturated, while Ed becomes saturated. Hence, the influx J will at some point become 
sensitive to X, while the deactivation rate 7 will no longer be proportional to [X*]. However, 
in the zero-order regime considered here, the total substrate concentration [S]t is large as 
compared to the enzyme concentrations and the Michaelis-Menten constants. This means 
that as [E a ]x is raised, such that [X](0) decreases and [X*](0) increases, initially E a will 
remain fully saturated, while Ed will become saturated. This implies that there is a range 
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of E a concentrations, where the influx J is still constant, but the deactivation rate 7 is no 
longer proportional to [X*]. In this regime, the concentration of X* increases with decreasing 
diffusion constant. Indeed, in this range, where [EJt < [Ed]T, the spatially non-uniform 
network will respond stronger than the spatially uniform network (see Figs. [3] and Eh-c). 

The significance of the diffusion term for [E^X*] in Eq. impedes a transparent analytical 
derivation of [X*](x) in this regime (see Text SI). However, under the condition that the 
influx J is constant, we can prove that the total amount of X* must decrease with increasing 
diffusion constant when 7 is no longer proportional to [X*]. The proof can be found in Text 
SI. 

Here we give a more intuitive explanation for the observation that a network in which 
the enzymes are spatially separated can respond stronger than a network in which the 
enzymes are similarly distributed in space. Ultimately, it is a consequence of the non-linear 
enzyme-substrate binding curve and the resulting hyperbolic dependence of the deactivation 
rate 7 on [X*] (see Eq. IT7l) . More specifically, this effect can arise when the diffusion 
constant is low and/or the deactivating enzyme operates close to, but not at, saturation 
in the uniform system; in this uniform system, X* is distributed evenly through the cell 
and all particles X* experience the same deactivation rate \i. In the spatially non-uniform 
system, [X*] is higher near the pole. If all the deactivating enzyme molecules would operate 
in the linear regime, i.e. if all deactivating enzyme molecules would not be saturated, then 
all particles X* would still experience the same degradation rate \i\ in this scenario, the 
increase in the number of X* particles close to the pole would precisely balance the decrease 
in the number of X* particles further away from the pole, as compared to the uniform 
network. However, if the concentration of the deactivating enzyme with respect to that of 
its substrate is lower, i.e. if the enzyme operates close to saturation, then the scenario can 
arise that the deactivating enzyme molecules near the pole become saturated (Fig. b), 
while in the corresponding uniform network they are not. In this scenario, the X* particles 
that are located close to the pole in the non-uniform network experience a lower effective 
deactivation rate than the X* particles in the spatially uniform network. This will enhance 
the response of the non-uniform system as compared to that of the uniform system. 

Strong activation 

We now discuss the effect of the diffusion speed on the concentration profiles of X* when 
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[E a ]T > [EJt (see Figs. [5]d-f ) . In this regime, [X] is low and [X*] is high. This reverses 
the saturation behavior of the antagonistic enzymes: while in the weak-activation limit E a 
is saturated and Ej is unsaturated, now E a is unsaturated and E^ is fully saturated. This 
also reverses the sensitivities of the influx J and the deactivation rate 7 to changes in the 
substrate concentration. Indeed, in the strong-activation regime not the influx J is constant, 
but rather the deactivation rate: 7 = /^[E^t (see Eq. H71) . 

In the limit that the deactivation rate 7 is constant, Eq. [15] can be solved in steady state 
(see Text SI). The solution is 

[X*](x) =co + Cl x + ^C2X 2 , (19) 

where c 2 = fc 6 [E d ] T /A c x = -k 6 [E d } T L/D, and c = [X*](0) = [S] T - [E d ] T (l + k 6 /k 3 (l + 
i^M,a/[E a ]T))- If is seen that in the high-activation regime, the concentration profile decays 
algebraically, rather than exponentially, as in the limit of weak activation. This is precisely 
because in the high-activation regime the total deactivation rate 7 is constant in space, while 
in the weak-activation limit 7 is proportional to the concentration of X*, which varies in 
space. In fact, in the weak-activation limit the deactivation rate per particle is constant in 
space and equal to \i. In contrast, in the strong-activation regime, the deactivation rate 
per particle is not only lower than fi on average, but also varies in space: the higher [X*] 
as compared to [E<i]t (which is constant in space and sets the total deactivation rate), the 
lower the deactivation rate per particle; activated particles close to the pole thus experience 
a lower deactivation rate and hence travel further on average before they are deactivated. 

The expression for cq = [X*](0) reveals that as [EJt increases, the concentration of X* 
close to the pole where E a is located, increases. In the limit that [E a ] — ► 00, [X*](0) — > [S]t — 
2[E d ]x ~ [S]t, which means that close to the pole of the cell where the activating enzyme 
is located, all the substrate X is converted into X* and EjX* (see Fig. [5H). Importantly, 
Fig. [5U also shows that as the distance from the pole increases, the fraction [X*] t (x)/[S]t(x) 
decreases, even in the limit that [E a ] T >> [E d ] T . When [E a ] T [E d ] T , all the substrate 
molecules at the pole will indeed be modified. However, these molecules will then diffuse 
away from the pole into the cytoplasm, where they can be demodified by the deactivating 
enzyme molecules, but not remodified. Hence, when the activating enzyme is spatially 
separated from the deactivating enzyme, it will never be possible to convert all the substrate 
molecules in the system (see Fig. [2]). This is in marked contrast with the situation in 
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which the activating and deactivating enzymes are not spatially separated. In this case, all 
substrate molecules can be converted into X* when [EJt ^> [Ed]T (see Fig. [2]). 

The expression for cq = [X*](0) also reveals that in the limit that [EJt 3> [EJtj the 
concentration of X* at x = is independent of the diffusion constant. However, while 
[X*](0) does not depend on the diffusion constant, the rate at which [X*](x) decays with the 
distance from the pole, does depend on it. Eq. [291 with c\ = —ke\Ed]TL/D, shows that the 
concentration profile of X* decays more slowly when the diffusion constant increases (see also 
Fig. [5H). These two observations, when taken together, imply that the total concentration 
of X* in the whole system increases with increasing diffusion constant. This can be verified 
by integrating Eq. |2"91 over the length of the cell, which gives [X*] ce11 ~ a — b/D, where a and 
b are positive constants. 

These results can be understood by comparing the influx of X* with the efflux of X*. 
When [EJt > [E(j]t, the deactivation rate is constant and hence independent of the diffusion 
constant. Since the total deactivation rate of X* is independent of the diffusion constant, 
the total influx of X*, which in steady state must balance the total efflux by deactivation, 
is also independent of the diffusion constant. The influx of X* depends on [E a X] and thus 
on the concentration of X* at x = 0, as discussed above. Hence, the concentration of X* 
at x = must be independent of the diffusion constant. A more intuitive explanation is as 
follows: As the diffusion constant increases, the X* molecules will diffuse away from the pole 
more rapidly. This would tend to lower the concentration of X* at x = 0. However, this 
process is accompanied by an increase in the flux of X towards the pole ([S]tO*0 is constant); 
because in the strong activation limit E a is unsaturated, this would tend to increase [E a X] 
and thereby the influx of X*, which would raise the concentration of X*. In steady state, 
these processes balance each other such that the concentration of X* at contact does not 
depend on the diffusion constant. However, while [X*](0) does not change with the diffusion 
constant, the X* molecules do diffuse away from the pole more rapidly when the diffusion 
constant increases. This means that the total concentration profile of X* must increase 
with increasing diffusion constant (see Fig. [3]). Indeed, only in the limit that D —> oo and 
[EJt ^> [Ed]T ; can all the substrate molecules be converted into X* (see Fig. [3]and Fig. [5H). 
In the strong-activation limit, spatially separating the antagonistic enzymes thus always 
weakens the response, in contrast to the behavior in the weak-activation limit. 

It should also be noted that the decay length of the concentration profile of X*, given by 
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c\ = — ke[E a ]TL/ D, does not only depend upon the diffusion constant, but also upon the 
activity of the deactivating enzyme. In the spatially non-uniform system, the maximum 
response (i.e. the response when [E a ] T ^> [E d ] T ) decreases as the catalytic activity of 
the deactivating enzyme, k$, increases. The reason is that the X* molecules will travel 
a shorter distance before they are deactivated, when the deactivation rate is higher. 
The extent to which spatially separating the enzymes weakens the maximum response 
thus depends upon both the diffusion constant and the deactivation rate of the X* molecules. 



Space-dependent amplification 

Fig. shows that if the activating enzyme is localized at one pole of the cell, while the 
deactivating enzyme can freely diffuse through the cytoplasm, the response of the network 
will depend upon the position in the cell. As can be deduced from Figs. H] and [5j [EdX*] 
depends significantly on the position in the cell when [E a ] T < [E d ] T . When [E a ] T > [E d ] T , 
however, [E d X*] becomes virtually independent of the position x, because then all the de- 
activating enzyme molecules are saturated. The opposite trend is observed for [X*]: when 
[E a ] T < [E d ] T , [X*] is low everywhere in the cell, while if [E a ] T > [E d ] T , [X*] strongly depends 
upon the position in the cell. The reason is, as discussed in the previous section, that even 
when [EJt ^> [E d ]x, not all X can be converted into X* if the two antagonistic enzymes are 
spatially separated. 

Interestingly, the average response of [E d X*] in the spatially non-uniform system is very 
similar to that in the system in which the two enzymes are not spatially separated. Yet, the 
response of [X*] does differ markedly between the two systems. This is a result of the strong 
non-linearity in the amplification mechanism of zero-order ultrasensitivity: because the ac- 
tivation and deactivation reactions are zero-order in the substrate concentrations [X] and 
[X*], respectively, even when fc 3 [E a X] is only marginally larger than fc 6 [E d X*], predominantly 
all X molecules will be converted into X* 4[]. 

Lastly, Fig. [6k shows that the inflection point of the dose-response curve depends on the 
position x in the cell. The inflection point shifts to higher [E a ]x/[E d ]x as the distance from 
the anterior pole increases; this effect becomes more pronounced as D decreases (data not 
shown). The fact that the inflection point depends on the position x is one of the principal 
reasons why the response in the spatially non-uniform system is weaker than that of the 
uniform system. 
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Push-pull networks in the linear regime 



Push-pull networks in living cells are not always in the zero-order regime 4|, |23| . In the linear 
regime, push-pull networks do not amplify signals, but can enhance the reliability of cell 
signaling by making it robust against fluctuations in the concentrations of the components 



due to noise in gene expression [24|]. It is therefore meaningful to study how the input-output 
relation of a push-pull network in the linear regime depends upon the spatial distribution of 
the antagonistic enzymes. A push-pull network in the linear regime is given by: 

E a + X^E a + X* (20) 
E d + X*^E d + X (21) 

The steady-state concentration profiles for these linear push-pull networks can be derived 
analytically. 

The principal result is that for push-pull networks that are in the linear regime, spatially 
separating the antagonistic enzymes always weakens the response. This can be seen by 
comparing the response curve for [S]t = OAKm in Fig. [2^ with that in Fig. [2b. The 
reason why for linear networks spatially separating the enzymes reduces the response 
in the strong-activation limit is the same as that for zero-order networks. The reason 
that, in contrast to zero-order networks, also in the weak- activation limit the response is 
weakened, is more subtle. In zero-order networks that are in the weak-activation limit, E a 
is saturated and, consequently, the influx J is independent of the concentration of X at 
the pole. In linear networks, E a is unsaturated and the influx J is proportional to [X](0). 
As D decreases, [X*](0) tends to increase and [X](0) tends to decrease ([S]t(#) is constant 
in space). Because in the linear regime J is proportional to [X](0), this would lower the 
influx of X*, which, in turn, would lower the concentration of X*. Spatially separating the 
antagonistic enzymes thus amplifies weak signals if the push-pull network operates in the 
zero-order regime, but not in the linear regime. 



Discussion 

In a push-pull network that operates deeply in the zero-order regime, the activation rate 
is given by A; 3 [E a ] T , while the deactivation rate is given by /c 6 [E d ] T ; both rates are thus 
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independent of the substrate concentration. If both enzymes are uniformly distributed, 
or colocalized, then essentially all substrate molecules will be activated when /^[EJt > 
fc 6 [E d ] T , while they will be predominantly deactivated when fc 3 [E a ] T < fc 6 [E d ] T . To drive the 
modification reactions to completion, it is indeed essential that the antagonistic enzymes 
are not spatially separated. If the antagonistic enzymes are separated, then the enzyme 
with the lower global activity, can locally still have a higher activity than the other enzyme. 
More in general, spatially separating the enzymes means that the balance between activation 
and deactivation depends upon the position in the cell, and this "smearing" of the response 
always tends to reduce the sharpness of the global response curve. 

If information about changes in the environment has to be transmitted, then the gain — 
the change in the output divided by the change in the input — is a critical quantity. In fact, 
the maximum gain is then usually the most relevant quantity, because signaling networks 
are often tuned to this point of maximum gain: the input-output function of a module 
and the concentration of its input signal are often optimized with respect to each other. 
The intracellular chemotaxis network of E. coli provides a clear example: the steady-state 
intracellular concentration of the messenger CheYp is around 3/xM, which is precisely the 
concentration at which the flagellar motors respond most strongly. Our analysis shows that 
from the perspective of signal amplification, the best strategy is to either colocalize the 
antagonistic enzymes or to uniformly distribute them in space: spatially separating the 
enzymes always weakens the maximum response. 

Nevertheless, as mentioned in the introduction, spatial gradients of messenger proteins are 
often observed. Indeed, maximizing the gain is not the only design principle in cell signaling. 
Firstly, while in some cases, such as E. coli chemotaxis, the signal has to be transmitted to 



a large number of places throughout the cell's cytoplasm or membrane [12], in other cases 
the signal has to be transmitted to distinct regions, such as the nucleus, or be confined to 
a small region near the membrane, as in the yeast pheromone response where the shmoo 
tip has to be formed locally; in this scenario, spatial gradients might be important, since 
they allow the cell to confine signaling to a narrow domain below the cell membrane p, [jij] . 
Secondly, a sharp response may not always be desirable. In order to respond strongly to 
changes in the input signal over a broad range of input signal strengths, the cell does not only 
need a sharp response curve, but it also needs to develop elaborate adaptation mechanisms 
that can reset the network to the point of maximum gain. In E. coli, for instance, the 
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methylation and demethylation enzymes CheR and CheB continually adjust the activity 
of the receptor cluster, such that the steady-state intracellular CheYp concentration is at 
3/iM. A weaker response curve, however, would allow the cell to have a reasonable working 
range without adaptation mechanisms. In this scenario not only the maximum gain would 
be important, but, in fact, the full response curve. Thirdly, it might not always be possible 
to maximize signal amplification by optimizing the input-output function of a module with 
respect to its incoming signal, because, for instance, the downstream module also has to 
respond to other incoming signals, while the signal also has to act on other downstream 
modules; the yeast MAPK pathways, which exhibit cross talk, provide a prominent example 
of such a scenario. It seems likely that in this case the full response curve, with the absolute 
concentrations of the components, is important. In this context it is interesting to note 
that spatially separating the antagonistic enzymes weakens strong signals by reducing the 
maximum output signal (Figs. [3]and[5ti), while it can enhance weak signals if the network 
operates in the zero-order regime (Figs. [3] and Eld) . This dependence of the input-output 
relation on the spatial distribution of the antagonistic enzymes could be exploited by cells 
to relay different environmental signals specifically. 

The analysis performed here is essentially a mean-field analysis. It is assumed that the 
concentrations are large and that fluctuations can be neglected. However, in the living cell, 
the concentrations are often low, which means that fluctuations can be important. This is 
particularly relevant for push-pull networks. Their high gain not only amplifies the mean of 
the input signal, but will also amplify the noise in the input signal [si, i25 . 
when the modification reactions become more zero order, the intrinsic fluctuations of the 
push-pull network, i.e. noise resulting from the modification reactions themselves, will also 

n 

increase |5[. In fact, it has been shown that when push-pull networks operate deeply in 
the zero-order regime, fluctuations can lead to a bimodal response Q. All these analyses 
of the effect of noise on the amplification mechanism of zero-order ultrasensitivity have 
been performed under the assumption that the enzymes are uniformly distributed in the 
cytoplasm. It would clearly be of interest to study the effect of enzyme (co)-localization on 
the noise characteristics of push-pull networks. 

Finally, could our predictions be tested experimentally? To test our predictions, one 
would ideally like to perform an experiment on a system with a canonical push-pull net- 
work in which all the parameters — concentrations of components, rate constants, diffusion 



26| . Moreover, 
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constants — are kept constant, except for the spatial location of one of the enzymes. This 
clearly seems a very difficult experiment to perform, and to our knowledge, no such experi- 
ment has been performed yet, with the possible exception of the experiment by Vaknin and 
Berg 1(3]. Vaknin and Berg studied the effect of phosphatase localization on the response 
of the intracellular chemotaxis network of E. coli cells. This network has a topology that is 
very similar to that of the canonical push-pull networks considered here, and it is believed 
that in the wild-type cells both the kinase and the phosphatase are localized at the cell pole. 
Vaknin and Berg compared the response of wild-type cells to that of mutant cells, in which 
the phosphatase was mutated such that it freely diffuses in the cytoplasm. They found that 
the spatial distribution of the phosphatase can have a marked effect on the the sharpness 
of the response, which seems to support the principal conclusion of our analysis. We would 
like to emphasize, however, that to assess the importance of the spatial distribution of the 
antagonistic enzymes in a push-pull network, a careful, quantitative analysis of the network 
is required. First of all, our analysis shows that both the quantitative and qualitative conse- 
quences of enzyme localization, depend upon the regime in which the network operates. For 
instance, our calculations reveal that if the activation rate is independent of the messenger 
concentration, and if the deactivation rate is linear in the messenger concentration, then the 
localization of the phosphatase should have no effect at all on the response curve. Secondly, 
it is quite possible that in the mutant cells not only the spatial distribution of the enzymes 
is different, but also their expression level, or even other parameters such as rate constants. 
In fact, experiments by Wang and Matsumura suggest that the activity of the phosphatase 
in the E. coli chemotaxis network is enhanced at the receptor cluster [27]. Clearly, differ- 
ent rate constants would also tend to change the response curve of the mutant cells with 
respect to that of the wild-type cells. To elucidate the effect of enzyme localization on the 
dose-response curve of a network, thus requires quantitative experiments and quantitative 
modeling. In a future publication, we will present a detailed analysis on the importance of 
phosphatase localization in the chemotaxis network oiE. coli. 

We thank Howard Berg, Dennis Bray, Viktor Sourjik and Ady Vaknin for useful dis- 
cussions and Rhoda Hawkins, Ady Vaknin and Viktor Sourjik for a critical reading of the 
manuscript. This work is part of the research program of the "Stichting voor Fundamenteel 
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LIST OF FIGURES 



1. A push-pull network. Two enzymes, E a and Ed, covalently (de)modify the components 
X and X*, respectively. The activating enzyme E a provides the input signal, the 
unmodified component X is the detection component and the modified component X* 
provides the output signal. 

2. The input-output relation of the push-pull network shown in Fig. [T] as a function of the 
total substrate concentration [S]t, for the case in which all components are uniformly 
distributed in space (a) and for the case in which the activating enzyme is located 
at one end of the cell, while the other components can diffuse freely through the cell 
(b). Here, [X*] T /[S] T = dx[X*} T (x) / dx[S} T (x). In (a) and (b), [E d ] T = 0.5 fM, 
Ky[,a = Km,6 = 0.5/iM, and k$ = k§ = 25s -1 . In (b), the diffusion constant is 
D = 10/im 2 s _1 . The inset shows the logarithmic gain g = din [X*]x/91n [E a ] T . It is 
seen that the sharpness of the response increases markedly with increasing substrate 
concentration when all the components are uniformly distributed in space (a), but 
much less so when the activating enzyme E a is located at one pole of the cell, while 
the deactivating enzyme Ed is distributed in the cytoplasm. When both enzymes E a 
and E<j are located at one pole, the steady-state dose-response curve is identical to 
that in (a). 

3. The input-output relation of a network in which the activating enzyme is located at 
one pole, while the other components can freely diffuse in the cytoplasm, for different 
values of the diffusion constant D (in /im 2 s _1 ) of the cytoplasmic components. The 
inset shows the logarithmic gain g = din [X*]x/91n [EJt- It is seen that the gain of the 
push-pull network strongly increases with increasing diffusion constant. If D — > oo, the 
dose-response curve approaches that of the push-pull network in which the components 
are uniformly distributed in space (and that of the network in which the enzymes are 
colocalized). The total substrate concentration is [S]t = 20/iM, the total concentration 
of the deactivating enzyme is [EJt = 0.5/iM, the Michaelis-Menten constants are 
Km,& = ^M,d = 0.5/iM, and the catalytic rate constants are k% = k§ = 25s" 1 . 

4. The concentration profiles for X* (a) and EdX* (b) in a push-pull network in which 
the activating enzyme is located at one pole of the cell, while the other components 
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are distributed in the cytoplasm, for three different concentrations of the activating 
enzyme. For all curves, [S] T = 20/iM, [E d ] T = 0.5/dVl, Km,& = K MA = 0.5 fM, 
h = k 6 = 25s _1 , and D = 10/im 2 s~ 1 . 

5. Profiles of [X*] (a and d), [EdX*] (b and e) and [Ed] (c and f). a - c: Low concentration 
of activating enzyme, [EJt = 0.5 [E^t; d-f: High activating enzyme concentration, 
[EJt = 1.5 [Ed] t- For the other parameter values, see Fig. HI 

6. Dose-response curves of the push-pull network in which the activating enzyme is lo- 
calized at one pole of the cell, while the other components diffuse in the cytoplasm, 
for different positions in the cell (x = corresponds to the black left most curve, 
while x = 3/mi corresponds to the black right most curve); a) profiles of [X*] and 
b) profiles of [EdX*]; note that the response becomes sharper further away from the 
pole. The green curves correspond to the average or integrated response of the non- 
uniform system, while the red curves correspond to the uniform system. The inset 
shows the logarithmic gain g = din [X*]/<91n [EJt at the respective positions in the 
cell (x = 0, 1, 2, 3, /mi). For the parameter values, see Fig. HI 
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Text SI 



The spatio-temporal evolution of the concentration profiles of a push-pull network where 
the activating enzyme is located at one pole of the cell, while the deactivating enzyme 
freely diffuses through the cytoplasm, is given by Eqs. 9-14 of the main text. To derive the 
spatio-temporal evolution of [X*] as given by Eqs. 15-17 of the main text, we have to make 
the assumptions that the formation of the enzyme-substrate complexes is fast and that the 
diffusion term for [E^X*] in Eq. 14 can be neglected. The steady-state solutions of Eqs. 12 
and 14 can then be combined with Eq. 9 to yield Eq. 15 of the main text. We now discuss 
the steady-state solution of Eq. 15 in the limits of weak-activation and strong-activation 
separately. 



Weak activation 

In the limit that [E a ]x <C [E^t? the diffusion term in Eq. 14 of the main text will be small 
compared to the other terms in steady-state. Combining the steady-state solution of that 
equation with those of Eqs. 9 and 12 of the main text, yields the following equation for X* 
in steady state: 

In the limit that [EJt <C [Ed]T) [X] will be large and [X*] will be small. In this limit, the 
above equation reduces to 

~ -h[W(x) + k 6 /K M4 [E d } T [X*}. (24) 

This equation can be solved with the boundary conditions 

d[X*](x 



D d[X*](x) 



dx 



J. 



dx 



= 0, (25) 

x=L 

where J = /^[EJt and L is the length of the cell. Defining the effective deactivation rate 
of X* to be [i = /c6/-^M,d[Ed]T and the characteristic decay length of X* to be A = y/D/ fi, 
then, if L ^> A, the solution of the above equation is 

[X*](*) = ^exp(-*/A). (26) 
30 



This equation predicts that the total amount of X* in the whole cell is independent of 
the diffusion constant D. This equation holds only in the limit that the diffusion term for 
[EdX*] in Eq. 14 of the main text can be neglected and if the deactivating enzyme is 
unsaturated, meaning that the total deactivation rate 7 is proportional to the concentration 
of X*. These assumptions are only accurate when [EJt <C [EJt- 

As [EJt increases, [X*] increases, and close to the pole where the activating enzyme re- 
sides, the deactivating enzyme does become saturated with X*. This lowers the deactivation 
rate per particle, which in turn increases the concentration of X* with respect to that in a 
system where both enzymes are uniformly distributed in space. 

The significance of the diffusion term for [E^X*] impedes an analytical derivation of 
[X*](x). However, when [E a ] T < [E d ] T and [X] is large, [E a X] is essentially constant 
and independent of the diffusion constant. Under this condition we can prove that 
the total amount of X* must decrease with increasing diffusion constant, implying that 
in this regime the uniform network will actually respond weaker than the spatially 
non-uniform network. Integrating Eq. 14 over the length of the cell reveals that in 
steady state k± dx[Ed](x)[X*](x) = (k^ + k§) dx\E^X*}{x). Combining Eq. 9 and 
Eq. 14 reveals that in steady state /cs[E a X] = k 6 dx[E d X*](x). Hence, in steady 
state [E a X] = c dx[Ed](x)[X*](x), where c is a constant. This means that if [E a X] is 
independent of the diffusion constant, also dx[E d ](x)[X*](x) must be independent of 
the diffusion constant. Now, when the diffusion constant decreases, the profiles [Ed](x) 
and [X*](x) vary more strongly in space and become less overlapping, as Figs. 5a and 
5c show. The decrease in overlap between [E a ](a;) and [X*](x) would tend to decrease 
J Q L c/x[Ed](x)[X*](x). To compensate for this, the total amount of X* must increase with 
decreasing diffusion constant, when [E a ]x < [E<i]t an d [E a X] is constant. 



Strong activation 

Combining Eq. 9 with Eq. 14 yields, in steady state: 

D^l = - + fc 6 [E d X*](x). (27) 

This equation for the steady-state concentration profile of X* is exact. Nevertheless, the 
diffusion term for [E^X*] impedes a transparant analytical solution. However, if [EJt > 
[E d ] T and D is not too low, then essentially all of the deactivating enzyme is saturated, and 
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[EdX*](x) = [E<i]t(^) = [Ed](a^) + [EdX*] (a;) (see Fig.5e of main text). Moreover, combining 
Eqs. 13 and 14 of the main text reveals that [E(i]t(^) is constant in space, if, as assumed here, 
the diffusion constants of the enzyme Ed, and that of the enzyme bound to its substrate, 
EdX*, are the same. Hence, in the limit that [E a ] is high and Ed is saturated, the second 
term on the right-hand-side of the above equation is zero, the third term is constant, and 
the equation reduces to 

D [ d J 2 { } = -k 3 [E a X]5(x) + k 6 [E d ] T . (28) 

This equation can be solved with the boundary conditions in Eq. [25] with J = fc 3 [E a X]. The 
solution is 

[X*](x) = c + Cl x + ^c 2 x 2 , (29) 

where c 2 = ke[E^/D, obtained from the solution of Eq. [28] in the domain < x < L, 
and c\ = — kelE^L/D, obtained from the boundary condition at x — L. The coefficient 
Co = [X*](0) can be obtained from the boundary condition at x = 0: /ce[Ed]T-^ = fca[E a X]. 
First, we note that the total substrate concentration in the cytoplasm, [S] c (x) = [X](x) + 
[X*](x) + [EdX*](x) = [S] c , is constant in space, because the diffusion constants of all the 
diffusing components are equal. Hence, the total substrate concentration in the whole cell 
is [S] T = [S] c + [E a X]/L = [S] T = [X](x) + [E a X]/L + [X*](x) + [E d ] T (where we have used 
that [EdX*] (a;) ~ [Ed]x)- The concentration [X](0) can be obtained by solving Eq. 12 of the 
main text in steady state, which gives [X](0) = i^M,a[E a X]/(L[E a ]T — [E a X]). Combining 
these expressions with the boundary condition at x — 0, and using that when [EJt is high, 
[E a ] T » [EaX], yields c = [X*](0) = [S] T - [E d ] T (l + h/k 3 (l + ir M , a /[E a ] T )). 
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